A new parallel strategy for two-dimensional incompressible flow simulations using pseudo-spectral methods

摘要

A novel parallel technique for Fourier-Galerkin pseudo-spectral methods withapplications to two-dimensional Navier-Stokes equations and inviscid Boussinesqapproximation equations is presented. It takes the advantage of the programmingstructure of the phase-shift de-aliased scheme for pseudo-spectral codes, andcombines the task-distribution strategy [Yin, Clercx and Montgomery, Comput.Fluids, 33, 509 (2004)] and parallelized Fast Fourier Transform scheme. Theperformances of the resulting MPI Fortran90 codes with the new procedure on SGI3800 are reported. For fixed resolution of the same problem, the peak speed ofthe new scheme can be twice as fast as the old parallel methods. Theparallelized codes are used to solve some challenging numerical problemsgoverned by the Navier-Stokes equations and the Boussinesq equations. Twointeresting physical problems, namely, the double-valued $\omega$-$\psi $structure in two-dimensional decaying turbulence and the collapse of the bubblecap in the Boussinesq simulation, are solved by using the proposed parallelalgorithms.